Incommensurate spin density modulation in a copper-oxide chain compound 

with commensurate charge order 
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Neutron diffraction has been used to determine the magnetic structure of NasCusOio, a stoichiometric com- 
pound containing chains based on edge-sharing Cu04 plaquettes. The chains are doped with 2/5 hole per Cu 
site and exhibit long-range commensurate charge order with an onset well above room temperature. Below 
Tjv = 23 K, the neutron data indicate long-range collinear magnetic order with a spin density modulation 
whose propagation vector is commensurate along and incommensurate perpendicular to the chains. Competing 
interchain exchange interactions are discussed as a possible origin of the incommensurate magnetic order. 



PACS numbers: 75.40.Gb, 63.20.-e, 78.30. -j, 75.50.-y 

The interplay between the magnetic and electric properties 
of copper oxides has recently been the subject of intense re- 
search activity. For instance, states with collinear and non- 
collinear magnetic order are currently under discussion in 
the contexts of ferroelectricity in undoped copper-oxide chain 
compounds [1] and of the anomalous transport properties of 
underdoped high temperature superconductors [2, 3], Since 
theoretical methods are well established in one dimension 
(ID), compounds with quasi- ID electronic structure are par- 
ticularly suitable as model systems to obtain a detailed under- 
standing of this interplay. However, research on the effect of 
doping in copper-oxide chain compounds has been limited by 
the scarcity of materials that support a significant density of 
holes on the chains. Most of the attention has been focused on 
the "telephone number compounds" (La,Sr) 1 4_ 2 .Ca ;c Cu2404i 
(LSCCO), which contain chain and ladder systems based on 
edge-sharing CUO4 square plaquettes [4, 5, 6, 7, 8, 9, 10]. 
Experiments have revealed intricate charge and spin ordering 
patterns on the chain subsystem, which depend strongly on 
the hole content. However, complications originating from the 
presence of two distinct electronically active subsystems with 
different hole concentrations and from the random potential of 
substituents (for x ^ 0) partially mask the intrinsic behavior 
of the copper-oxide chains. Moreover, recent work has shown 
that the magnetic properties of this material are strongly in- 
fluenced by an incommensurate structural modulation arising 
from a mismatch of different units constituting the crystal lat- 
tice [8, 9, 10]. Ca 2+l; Y2_ a; Cu50io (CYCO), a class of materi- 
als containing only copper-oxide chains, also exhibits a com- 
plicated structural modulation unrelated to charge ordering 
[11, 12]. In addition, the magnetic properties of doped chains 
in this material appear to be influenced to a large extent by 
substitutional disorder and/or oxygen non-stoichiometry [13]. 

Naa;Cu02, a recently synthesized family of compounds 
with very low chemical disorder, offers new perspectives in 
this regard [14, 15, 16]. This material consists entirely of 
electronically inert Na + ions and chains built of edge-sharing 
CUO4 plaquettes similar to those in LSCCO and CYCO (Fig. 



1). Holes donated to the chains by the Na ions form long- 
range ordered superstructures that are generally incommensu- 
rate with the Na sublattice. However, in contrast to other cop- 
per oxides with dopable chains, incommensurate structural 
modulations due to purely steric constraints are not present, 
so that commensurate charge order can be established if x is 
a rational number. By carefully adjusting the chemical syn- 
thesis conditions, a state with x — 1.60, corresponding to a 
hole filling factor of 2/5 on the chains, has recently been stabi- 
lized [16]. The stoichiometric compound created in this way, 
NagCusOio, is a unique testing ground for theories of mag- 
netism in doped copper oxides, without complications aris- 
ing from substitutional disorder and/or incommensurate lat- 
tice distortions. 

We have determined the magnetic structure of NagCusOio 
below its Neel temperature T/v = 23 K by neutron powder 
diffraction. We find that the spins are collinear and exhibit an 
incommensurate spin density modulation that is unusual for 
magnetic insulators. A possible origin is a network of com- 
peting interchain exchange interactions. An investigation of 
the mechanisms stabilizing this state compared to the heli- 
coidal states found in the undoped analogues (Na,Li)Cu202 
[17, 18, 19] may provide important insights into the compe- 
tition between collinear and noncollinear magnetism in other 
copper oxides of topical interest [1,2, 3]. 

Powder samples of NagC^Oio with weight ~ 4.5 g were 
synthesized in a single batch as described previously [14]. 
Their magnetic susceptibility was found to be in good agree- 
ment with prior reports [14, 15]. In particular, a magnetic tran- 
sition temperature Tjv = 23 K was obtained by analyzing the 
derivative of the magnetization as a function of temperature. 
As NagCusOio is sensitive to air, the samples were sealed in 
air-tight vanadium cans, which were loaded into a helium flow 
cryostat. The neutron diffraction data were taken at the Labo- 
ratoire Leon Brillouin in Saclay, France. In order to determine 
the nuclear structure, we used the high-resolution diffractome- 
ter 3-T-2 with a neutron wavelength of 1.23 A. The data for 
the magnetic structure determination were taken on the high- 



2 




Cul Cu3 Cu2 

f • ♦ 




0Cu+2 
©Cu+3 



FIG. 1: (a) Nuclear unit cell and spin structure of NagCusOio- (b) 
Cut of the lattice structure along the ab plane, showing inequivalent 
copper sites and superexchange parameters, as explained in the text. 



flux cold-neutron diffractometer G-4- 1 with a neutron wave- 
length of 2.43 A. 

Fig. 2 shows a high-resolution powder diffraction pattern 
obtained at room temperature. The nuclear intensities for both 
samples were refined with the program Fullprof [20] on the 
basis of the monoclinic space group Cm. The unit cell param- 
eters a = 8.23492 ± 0.00014A, b = 13.92889 ± 0.00020A, 
c = 5.71324 ± 0.00010A, and (3 = 111.91 ± 0.004° and 
atomic positions obtained from the refinement agree with the 
results of earlier x-ray diffraction studies [14, 16]. The re- 
sulting diffraction pattern yields an excellent description of 
the experimental data (Fig. 2), as indicated by the goodness- 
of-fit-parameter Rp = 0.0854. This confirms that the sam- 
ple is chemically homogeneous, and that the lattice structure 
is commensurate. Another sample prepared under nominally 
identical conditions yielded substantially worse refinements, 
and data on a sample consisting of batches of powder mate- 
rial synthesized in different reactions could only be fitted by 
a superposition of different charge ordering patterns. These 
findings confirm that slight deviations from the ideal compo- 
sition result in incommensurately modulated structures [16]. 

The atomic positions within the unit cell are displayed in 
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FIG. 2: High-resolution neutron powder diffraction pattern of 
NagCusOio at room temperature. The black line shows the result 
of the Rietveld refinement discussed in the text. The positions of nu- 
clear Bragg reflections are indicated by green marks. The blue curve 
gives the difference between the calculated and measured intensities. 



Fig. 1. The unit cell comprises ten copper ions, which are 
organized in two parallel CuC>2 chains pointing along the b- 
axis. The chains are separated by Na ions. Four of the copper 
ions in each unit cell (Cu3 in Fig. lb) were found to exhibit 
bonding patterns characteristic of spinless Zhang-Rice singlet 
states with formal valence 3+ (Ref. 14). The Cu-O bond 
lengths of the remaining six copper ions indicate a valence 
state of Cu 2+ with spin 1/2. The Cu 2+ ions are located in two 
inequivalent sites, which are surrounded by two Cu 3+ ions 
(Cul), and one Cu 3+ ion and one Cu 2+ ion (Cu2), respec- 
tively. Nominally di- and trivalent copper ions are ordered 
in the sequence 2-2-3-2-3-2-2-3- ... along the chains (Fig. 1). 
The charge order is stable up to temperatures well above room 
temperature [14, 15]. 

Fig. 3 shows the low-angle segment of the high-flux pow- 
der pattern measured at T — 1.4 K taken to determine the 
magnetic structure. Two of the peaks shown vanish above 
Tyv = 23 K identified from the anomaly in the uniform sus- 
ceptibility and hence originate from magnetic scattering. The 
temperature dependence of the magnetic Bragg intensity is 
well described by a power-law fit without detectable round- 
ing near Tn (Fig. 4), which indicates homogeneous magnetic 
long-range order in the low-temperature phase. This conclu- 
sion is supported by the refinement described below, which 
shows that the widths of the magnetic Bragg reflections are 
limited by the instrumental resolution. 

The magnetic structure was refined using the software 
package Fullprof [20] based on the magnetic structure factor 
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FIG. 3: Low-angle segment of the high-flux powder diffraction pat- 
tern measured at 1.4K (red points). The black line shows the results 
of the Rietveld refinement discussed in the text. The upper (lower) 
green marks indicate the positions of nuclear (magnetic) Bragg re- 
flections. The blue curve gives the difference between the calculated 
and measured intensities. The inset shows the result of a refinement 
based on a commensurate spin structure. 



where h is the scattering vector, k is the propagation vec- 
tor of the magnetic structure, j enumerates the symmetry- 
inequivalent magnetic ions at positions fj in the magnetic unit 
cell, the index s runs over the magnetic symmetry operators, 
Mj S is an operator that transforms the Fourier components 
Sgj of the magnetic moments according to a given symme- 
try, {S|i} s generates the positions of all symmetry-equivalent 
magnetic ions in the unit cell, and ^£j S is a phase factor. The 

refinements are based on the isotropic form factor fj(h) of 
the Cu 2+ ion. The occupation factor Oj as well as the Debye- 
Waller factors, T- so and T JS , were set to unity. 

The positions of the magnetic Bragg peaks indicate an ap- 
proximate doubling of the unit cell along a and c, whereas 
the magnetic and nuclear unit cells coincide along the spin- 
chain axis b. A systematic shift away from scattering an- 
gles corresponding to commensurate Bragg reflections (inset 
in Fig. 3) reveals that the magnetic structure is incommen- 
surate. Note that this shift can be extracted with high confi- 
dence, because the reference lattice parameters are determined 
by many Bragg reflections over a wide range of scattering an- 
gles (Fig. 2). The propagation vector resulting from the re- 
finement is (-0.5 + a, 0, 0.5 - j3) with a = 0.089(3) and 
(3 = -0.030(1) (at 1.4 K). The Miller indices of the mag- 
netic Bragg reflections are shown in Fig. 3. The asymmetric 
lineshape of the magnetic peak at higher scattering angle is 
well explained as a consequence of the superposition of two 
resolution-limited, nearly coincident incommensurate Bragg 
reflections. 

Although only three inequivalent magnetic Bragg reflec- 
tions are visible, the diffraction pattern imposes strong con- 
straints on the magnetic structure. Because of the large 



FIG. 4: Integrated intensities of the (=p0.5=F«, 1, ±0.5±/3)- and (0.5- 
a, 0, 0.5+/3) magnetic Bragg reflections as a function of temperature. 
The line is the result of a power-law fit. 



unit cell and the incommensurate magnetic modulation, most 
possible spin arrangements generate additional Bragg reflec- 
tions with intensities well outside the experimental error bars, 
where none are observed. By far the best agreement with the 
data (Bragg R = 0.049, magnetic R — 0.114) was obtained 
based on a collinear spin structure in which the two Cu 2+ mo- 
ments on Cu2 sites directly adjacent along the chains are par- 
allel, whereas Cul and Cu2 moments separated by Cu3 ions 
are antiparallel. The magnetic moment on the Cu3 site was 
refined to zero, consistent with the Zhang-Rice singlet state 
inferred from the bond-length analysis of Ref. 14. The in- 
commensurate propagation vector perpendicular to the chains 
modulates the magnitude of the Cu 2+ moments. The modula- 
tion amplitude was refined as 0.84±0.10/zs on both Cul and 
Cu2 sites, consistent with a spin- 1/2 state. The spin direction 
resulting from the refinement is (0.86 ± 0.39, 0, 0.92 ± 0.07). 
An additional refinable parameter is the phase difference 
*i2 = 44.6 ± 4.5° of the modulation on Cul and Cu2 sites. 
The spin arrangement within the nuclear unit cell is displayed 
in Fig. 1 . The corresponding diffraction pattern is in excellent 
agreement with the data (Fig. 3). 

A comprehensive set of alternative collinear and non- 
collinear spin structures was also tested, but the resulting re- 
finements were unsatisfactory. In particular, the diffraction 
patterns of the circular helix structure that yields the best 
agreement with the data generates prominent Bragg reflec- 
tions (±0.5 ± a, 1, ±0.5 T P) at 22.9°, where the Bragg in- 
tensity vanishes within the error (arrow in Fig. 3). The corre- 
sponding magnetic i?-factor is 0.354, much worse than that of 
the collinear state. If the refinement is generalized to include 
elliptical helix structures, the length of the minor axis of the 
ellipse (0.061 ± 0.165/is along b) is consistent with zero, and 
the magnetic i?-factor does not improve significantly com- 
pared to the collinear state. Although a small noncollinear 
component due to helicity or canting cannot be ruled out, the 
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incommensurate modulation therefore predominantly affects 
the moment amplitude. 

The observed spin amplitude modulation is formally anal- 
ogous to spin density waves in metallic systems, but the in- 
sulating nature and robust charge order of Na 8 Cu 5 Oio im- 
ply that it cannot arise from a Fermi surface instability. We 
therefore discuss our data in terms of superexchange interac- 
tions between local magnetic moments, focusing first on the 
commensurate spin structure along the chain axis b. The spin 
alignment along this axis indicates a ferromagnetic nearest- 
neighbor (nn) exchange interaction Jp and an antiferromag- 
netic next-nearest-neighbor (nnn) interaction J» 2 (Fig. lb), 
in agreement with electronic structure calculations for edge- 
sharing copper-oxide chains [7, 15, 21, 23, 24, 25] and with 
the conclusions of experiments on a variety of undoped com- 
pounds including LiCu202 [17, 19] andNaCu2C>2 [18], which 
contain undoped chains with similar bond lengths and angles 
as the ones in Na 8 Cu 5 Oi . Since the Cu-O-Cu bond angle 
in the edge-sharing chain geometry is close to 90°, Jim is 
anomalously small, and the competing nnn coupling J|| 2 is 
comparable or larger in magnitude. The undoped spin sys- 
tems of (Li,Na)Cu2C>2 respond to the resulting frustration by 
forming incommensurate, helical magnetic order propagating 
along the chains [17, 18, 19]. In Na 8 Cu 5 Oio, charge order- 
ing lifts this frustration and gives rise to a commensurate spin 
structure along the chains, as predicted by model calculations 
[15]. 

The situation is different for interactions between different 
chains. For simplicity, we first ignore the small incommen- 
surability along c and consider the magnetic bonding pattern 
in the a6-plane (Fig. lb), including both interactions between 
directly adjacent chains within the same unit cell (Jj_i, J ±2) 
and interactions between nnn chains ( J13). The simplest ex- 
planation for the approximate doubling of the unit cell along a 
is that J13 is dominant and antiferromagnetic, leaving the in- 
teractions between nn chains frustrated. In principle, the spin 
system can respond to the frustration by establishing either 
noncollinear magnetic order, as observed in (Na,Li)Cu2C>2, 
or periodic spin-singlet correlations, as found in models of 
frustrated and/or doped ID [4, 6] and 2D [3] quantum anti- 
ferromagnets. An admixture of such correlations is a possible 
mechanism underlying the observed spin density modulation 
in NagCusOio- To obtain a crude estimate of the magnitude of 
the interchain interactions in the framework of this scenario, 
we consider collinear classical spins coupled by sinusoidally 
modulated exchange bonds with amplitudes shown in Fig. lb. 
Minimization of the exchange energy with respect to the phase 
shift ^12 of the modulation on the Cul and Cu2 sublattices 
then yields J±i/Jp = sin ^i 2 /2 sin(ifc x a + ^12) ~ —0.7, 
where k x is the component of the incommensurate propaga- 
tion vector along the a-axis determined by the competing in- 
teractions between nn and nnn chains. While a full quantum- 
mechanical calculation is required to assess the viability of 
this scenario, this simple estimate indicates that the exchange 



interactions between CUO2 chains along a are comparable to 
those within the chains, as observed for other copper-oxide 
chain compounds [22, 23, 24]. The smaller incommensura- 
bility along c suggests weaker exchange interactions in this 
direction. 

In view of the helicoidal states observed in (Li,Na)Cu 2 02, 
the collinear spin density modulation in Na 8 Cu 5 Oio may 
seem surprising. However, as other cuprates with undoped 
edge-sharing chains exhibit collinear spins [12, 22], the en- 
ergy balance between both types of order appears to be 
quite subtle. This is confirmed by ab-initio calculations 
[23, 25]. Anisotropic exchange [5, 22] and/or order-from- 
disorder mechanisms [26] may be responsible for tipping the 
balance towards collinear order in Na 8 Cu 5 Oio- 

We thank F. Bouree and B. Rieu for help during the mea- 
surements at LLB, P. Horsch, R. K. Kremer, and R Bourges 
for useful discussions, and E. Briicher for the susceptibility 
measurements. 



[1] See, e.g., S. Parke? al, Phys. Rev. Lett. 98, 057601 (2007); H.J. 

Xiang and M.-H. Whangbo, ibid. 99, 257203 (2007); Y. Naito 

et al, J. Phys. Soc. Jpn. 76, 023708 (2007). 
[2] See, e.g., A. Liischer, A.I. Milstein, and O.P Sushkov, Phys. 

Rev. Lett. 98, 037001 (2007). 
[3] See, e.g., S. Papanikolaou, K.S. Raman, and E. Fradkin, Phys. 

Rev. B 75, 094406 (2007); M. Vojta and S. Sachdev, Phys. Rev. 

Lett. 83, 3916 (1999). 
[4] M. Matsuda et al, Phys. Rev. B 59, 1060 (1999). 
[5] V. Kataev et al, Phys. Rev. Lett. 86, 2882 (2001). 
[6] R. Klingeler et al, Phys. Rev. B 73, 014426 (2006). 
[7] U. Schwingenschlogl and C. Schuster, Phys. Rev. Lett. 99, 

237206 (2007). 

[8] A. Gelle and M.B. Lepetit, Phys. Rev. Lett. 92, 236402 (2004). 
[9] M. v. Zimmermann et al, Phys. Rev. B 73, 115121 (2006). 
[10] A. Rusydi et al, Phys. Rev. Lett. 100, 036403 (2008). 
[11] PK. Davies, J. Sol. State Chem. 95, 365 (1991). 
[12] H. F. Fong et al., Phys. Rev. B 59, 6873 (1999); M. Matsuda, 

K. Ohyama, and M. Ohashi, J. Phys. Soc. Jpn. 68, 269 (1999). 
[13] M.D. Chabot and J.T. Markert, Phys. Rev. Lett. 86, 163 (2001); 

M. Matsuda et al, Phys. Rev. B 71, 104414 (2005); K. Kudo et 

al, ibid. 71, 104413 (2005). 
[14] M. Sofin et al, J. Sol. State Chem. 178, 3708 (2005). 
[15] P. Horsch etal, Phys. Rev. Lett. 94, 076403 (2005). 
[16] S. van Smaalen et al, Acta Cryst. B 63, 17 (2007). 
[17] AA. Gippius etal, Phys. Rev. B 70, 020406(R) (2004). 
[18] L. Capogna etal, Phys. Rev. B 71, 140402(R) (2005). 
[19] T. Masuda et al, Phys. Rev. B 72, 014405 (2005). 
[20] J. Rodriguez-Carvajal, PhysicaB 192, 55 (1993). 
[21] Y. Mizuno et al., Phys. Rev. B 57, 5326 (1998). 
[22] M. Boehm et al, Europhys. Lett. 43, 77 (1998). 
[23] H.J. Xiang et al, Phys. Rev. B 76, 22041 1(R) (2007). 
[24] S.L. Drechsler etal, Europhys. Lett. 73, 83 (2006). 
[25] S.L. Drechsler et al, Phys. Rev. Lett. 98, 077202 (2007). 
[26] Y.J. Kim et al, Phys. Rev. Lett. 83, 852 (1999). 



